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ABSTRACT 

We model the 1990 giant X-ray flare of the quiescent galaxy NGC 5905 as 
the tidal disruption of a star by a supermassive black hole. From the observed 
rapid decline of the luminosity, over a timescale of a few years, we argue that the 
flare was powered by the fallback of debris rather than subsequent accretion via a 
thin disk. The fallback model allows constraints to be set on the black hole mass 
and the mass of debris. The latter must be very much less than a solar mass to 
explain the very low luminosity of the flare. The observations can be explained 
either as the partial stripping of the outer layers of a low-mass main sequence star 
or as the disruption of a brown dwarf or a giant planet. We find that the X-ray 
emission in the flare must have originated within a small patch rather than over 
the entire torus of circularized material surrounding the black hole. We suggest 
that the patch corresponds to the "bright spot" where the stream of returning 
debris impacts the torus. Interestingly, although the peak luminosity of the flare 
was highly sub-Eddington, the peak flux from the bright spot was close to the 
Eddington limit. We speculate on the implications of this result for observations 
of other flare events. 

Subject headings: accretion, accretion disks — black hole physics — galaxies: 
nuclei 
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1. Introduction 

Tidal disruption of a star by a supermassive black hole in a galactic nucleus has been 
investigated by many authors. Early seminal investigations by Hills (1975); Lacy, Townes, 
& Hollenbach (1982); Rees (1988); Evans & Kochanek (1989) and Cannizzo, Lee, & Good- 
man (1990) were followed up with detailed analytical studies (Kochanek 1994; Khokhlov, & 
Melia 1996; Loeb & Ulmer 1997; Ulmer, Paczyhski, & Goodman 1998; Ulmer 1999) and nu- 
merical simulations (Khokhlov, Novikov, & Pethick 1993; Frolov et al. 1994; Marck, Lioure, 
& Bonazzola 1996; Diener et al. 1997; Kim, Park, & Lee 1999; Ayal, Livio, & Piran 2000; 
Ivanov & Novikov 2001). 

A tidal disruption event is expected to produce a luminous flare of electromagnetic 
radiation in the UV to X-ray band. Komossa (2001) argued that a convincing detection of 
a tidal disruption event must fulfill the following three criteria: (1) The event should be of 
short duration (a "flare"); (2) It should be very luminous (up to L max ~ 10 45 erg s _1 at 
maximum); (3) It should reside in a galaxy that is otherwise non- active (in order to exclude 
an upward fluctuation in the mass accretion rate of an active nucleus). 

R OSA T observations have led to the discovery of four X-ray flaring "normal" galaxies 
that fulfill the above criteria [summarized in Komossa (2001)]: NGC 5905, RXJ1242-1119, 
RXJ1624+7554, RXJ1420+5334; plus a possible fifth candidate: RXJ1331-3243. All these 
galaxies have exhibited large X-ray flares (peak luminosity up to ~ 10 44 erg s -1 ), corre- 
sponding to a significant amplitude of variability (up to a factor ~ 200), an ultra-soft X-ray 
spectrum (k B T hh 0.04 — 0.1 keV, where T bb is the temperature derived from a black-body 
fit to the spectral data), and no signs of Seyfert-like activity in the optical. These flares 
are promising candidates for tidal disruption events (Komossa 2001). A UV flare has also 
been detected at the center of a mildly active elliptical galaxy, NGC 4552, with a bolometric 
luminosity ~ 10 39 erg s" 1 . Since the luminosity is several orders of magnitude less than 
would be expected from complete tidal disruption of a star, the event has been explained as 
the tidal stripping of a stellar atmosphere (Renzini et al. 1995). 

Among the five galaxies listed above, long-term follow-up observations exist only for 
NGC 5905 (Komossa & Bade 1999). In this paper, we model the light curve of NGC 5905 
assuming that the flare was produced by the tidal disruption of a star by a supermassive 
black hole. We compare theoretical predictions with the observational data and constrain 
the properties of the progenitor star (e.g., its mass and radius) and estimate how much mass 
was accreted by the black hole. 
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2. The Light Curve of the Flare in NGC 5905 

A giant X-ray flare in NGC 5905 was detected by ROSAT'm July 1990 (Bade, Komossa, 
& Dahlem 1996). The observed peak X-ray luminosity exceeded 10 42 erg s _1 , and the X-ray 
spectrum was very soft (photon index T x ~ —4, or a very cool black-body with /c^Tbb « 0.06 
keV). Follow-up observations were made in 1990 and 1993 with the PSPC, and in 1996 with 
the HRI (Komossa & Bade 1999). The observational data are listed in Table 1. Assuming a 
Hubble constant H = 50 km s _1 Mpc -1 , the count rates may be converted to luminosities, 
which are plotted with error bars in Figure 1. The two downward arrows are upper limits. 
The light curve after the initial burst gradually declines with time. During the X-ray outburst 
as well as at later times, the optical brightness is constant to within observational errors 
(~ 0.2 m ). 

Fitting the late decline of the NGC 5905 light curve (the four last data points in Table 1) 
with a (t — £d)~ 5//3 law that is predicted by the fallback model of tidal disruption (see §3), 
and treating the largest luminosity observed at t — 1990.54 yr, viz. Lx = 4.47 x 10 42 ergs^ 1 , 
as a lower limit (since this was still in the rising part of the light curve), we obtain 

L x = (0.30 ±0.03) x lO^ergs -1 

where the time t is in units of years, and t D = 1990.36 ± 0.02 yr is the moment when tidal 
disruption happened. 

During the first epoch of observations at time t = 1990.54, the X-ray flux rose steadily 
(see Table 1), strongly suggesting that this epoch corresponded to the onset of the flare. For 
a start time of t 1 ~ 1990.54 yr, the total X-ray energy released during the entire outburst 
event is 

AE X = L x (t)dt « (4.5 ± 0.9) x 10 49 erg . (2) 

Jt! 

A black-body fit to the X-ray spectrum during the outburst gives a temperature A; B T bb « 
0.06 keV, and an intrinsic (0.1 — 2.4 keV) mean luminosity L x ~ 3 x 10 42 erg s^ 1 . Integrating 
the black-body spectrum over all energies (from zero to infinity), the bolometric luminosity 
is L R3 l.lLx- Thus, if the spectrum does indeed have a black-body form, the bolometric 
correction is not very significant. Then, the total mass associated with the released energy 
may be estimated to be 

AM = w w (2.5 ± 0.5) x 10" 4 M (^-) , (3) 



t- (1990.36 ± 0.02) yr 

1 y r 
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where e is the efficiency of converting mass to radiated energy. This mass estimate is similar 
to that for the UV flare in NGC 4552 (Renzini et al. 1995). Note that our estimate of AE X 
depends only weakly on the assumed power law index 5/3 in equation (1). This is because 
the observed light curve extends over a long enough time to capture most of the emission 
in the flare. The estimate of AM suffers from the uncertainty in the bolometric correction 
(which we have taken to be unity). 

In models of thin accretion disks around black holes, the color temperature Tbb of 
the emitted radiation is generally higher than the effective temperature T e fj, defined by flux 
F = o"T c g. The deviation of Tbb from T e g is caused by electron scattering and Comptonization 
(Ross, Fabian & Mineshige 1992; Shimura & Takahara 1993, 1995a, 1995b). Writing Tbb = 
f c T e s, the factor f c can be estimated from detailed self-consistent models of the radiative 
transfer in the disk. It is found that the value of f c depends on both the mass of the 
black hole (compare Shimura & Takahara 1993 and 1995b) and the mass accretion rate. 
Shimura & Takahara (1993) present results for a 1O 8 M black hole emitting at a radius of 5 
Schwarzschild radii, which is similar to the conditions we expect in our problem. Depending 
on the accretion rate, they obtain values of f c in the range 1.3 — 3.2 (see their Fig. 4 and 
the associated discussion), with the highest value of f c being obtained for near-Eddington 
rate accretion. In what follows, we retain the factor f c in the equations, and we substitute 
a value of 3 whenever we require a numerical estimate. 

Including the factor f c as defined above, we estimate the radius of the X-ray emitting 
region of the flare in NGC 5905 to be 



where we have used Lx ~ 3 x 10 erg s for the mean luminosity during the peak of 
the outburst (Komossa 2001). The radius we derive is surprisingly small; even after setting 
f c = 3, it is barely equal to the Schwarzschild radius of a 1O 7 M black hole, whereas an 
accretion disk is likely to have an effective radiating area 10 — 100 times larger. We return 
to this point in §5. 

3. Tidal Disruption of a Star by a Supermassive Black Hole 

The tidal disruption radius of a star in the vicinity of a Schwarzschild black hole is 




(4) 




(5) 
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where M H is the mass of the black hole, is the mass of the star, is the radius of 
the star, and /i is a dimensionless coefficient. The exact value of /x is uncertain, though 
physically we expect it to be of order unity. We assume /i — 1 in our modeling. 

The ratio of r T to r H , where r H = 2GM H /c 2 is the Schwarzschild radius of the black 
hole, is 



For a solar- type star with M* = M Q and = R & , we have tt < r# when Mh > 1.1 x 
10 8 M Q . In such a case the star will be swallowed whole by the black hole, without a tidal 
disruption event (Hills 1975). 

After a star is tidally disrupted, there are two distinct stages in the dynamical evolution 
of the debris, each with its own time development and luminosity: 

(1) The fallback stage (Rees 1988): A fraction of the material in the disrupted star 
remains gravitationally bound to the black hole and falls back to the pericenter, giving rise 
to a mass accretion rate (and a luminosity) evolving with time as ~ 

j-5/3 (Phinney 1989). 

The returning material shocks and then circularizes to form an orbiting torus around the 
black hole at a radius r c = 2rp (by angular momentum conservation), where rp is the 
pericentric radius. The total energy radiated during the fallback is given by GM H AM/2r c , 
where AM is the fallback mass. This energy can be less than the energy released in the 
subsequent viscous accretion stage discussed below. Nevertheless, since the fallback stage is 
relatively short-lived, it dominates the early luminosity of a disruption event. 

(2) The viscous accretion stage (Cannizzo, Lee, & Goodman 1990): The torus formed 
from returning debris gradually spreads inward and outward from r — r c under the action of 
viscosity, and gives rise to a mass accretion rate (and a luminosity) evolving with time approx- 
imately as ~ t~ 12 . The energy radiated during this stage is given by GMi{AM[(l/2r ms ) — 
(l/2r c )], where r ms is the radius of the last stable circular orbit (= 3 Schwarzschild radii for 
a non-spinning black hole). The accretion occurs on a viscous time scale, which is typically 
very long compared to the fallback time (Cannizo et al. 1990; Ulmer 1999; Appendix A of 
the present paper). Therefore, this stage is expected to last a long time, up to hundreds 
of years, depending on the magnitude of the viscosity in the disk, and the luminosity is 
significantly lower than that associated with the fallback stage (Appendix A). 

In the case of NGC 5905, the X-ray luminosity was seen to rise rapidly during the first 
epoch of observations at t — 1990.54 yr, indicating that the flare probably began around 
that time (see Table 1 and Fig. 1). The luminosity then dropped by more than a factor of 
100 between t = 1990.54 yr and t = 1996.89 yr. Such a rapid decline indicates that the 
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observed outburst of NGC 5905 could not have been due to the accretion phase, but must 
have corresponded to the fallback stage (if it indeed was a tidal disruption event). Therefore, 
in this paper we assume that the flare in NGC 5905 corresponds to the fallback stage of tidal 
disruption. 

Let us assume that the center of the star follows a nearly parabolic orbit (i.e., the binding 
energy of the star to the black hole is close to zero), and for definiteness let us assume that the 
pericenter of the orbit is at rp = ty, i.e. that the parameter rj defined below in equation (8) is 
equal to unity. Since the specific energy at the center of the star is E c = v 2 p /2 — GMh/tt = 0, 
the orbital velocity of the star at pericenter is vp — (2GMh/ttY^ 2 = c (r^/r^) 1 / 2 , and the 
spread in the specific energy of the disrupted debris, 2AE, is governed by the variation of 
the black hole gravitational potential across the star, and the spin-up of the star as a result 
of the tidal interaction (Rees 1988). If M H 3> M*, we have 



where k depends on the spin-up state of the star. If the star is spun up to the break-up 
spin angular velocity, we have k ~ 3. On the other hand, if the spin-up effect is negligible, 
then we have k ~ 1 (Rees 1988; Lacy, Townes, & Hollenbach 1982; Evans & Kochanek 1989; 
Ayal, Livio, & Piran 2000). Since the spin- up velocity is always C t/p, the energy used to 
spin up the star is much less than the orbital kinetic energy of the star, and the stellar orbit 
remains parabolic to a very good approximation. 

The spin-up of a star via tidal interaction is a complex process. In linear perturbation 
theory, the spin- up angular velocity is given by (Press & Teukolsky 1977; Alexander & Kumar 
2001; Alexander & Livio 2001) 



where 77 ~ 1 in our problem, uj p = vp/r P is the orbital angular velocity of the star at the 
pericenter, a is the stellar momentum of inertia in units of M*/? 2 , and T 2 is the second 
tidal coupling coefficient which is a function of rj. In equation (8) we have assumed that 
Mh ^> M*. For a main sequence star of mass 0.76M Q and radius O.75i? , Alexander & 
Kumar (2001) give a ~ 0.07 and T^l) ~ 0.06, which implies that uo s ~ 0.43cjp if 77 = 1. 
For an n = 1.5 polytrope star, they give a ~ 0.21 and T 2 (l) ~ 0.36, which corresponds to 
uj s 0.86u;p. Furthermore, they show that numerical simulations that include non-linear 
effects lead to a larger energy transfer from the orbit to the star and a larger spin-up than 
predicted by linear theory. For these reasons, we feel that it is likely a tidally disrupted star 
will be spun up to nearly the break-up limit, with oj s ~ u P and k ~ 3. We therefore assume 




(7) 




(8) 
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k = 3 for most of our analysis; however, for completeness we also present a few results for 
the non-spinning case k — 1. 

Part of the disrupted star has positive energy and escapes from the system, and part 
remains bound. The bound material follows a highly eccentric orbit and returns to the central 
black hole after completing one orbit. The flare begins when the most bound material, with 
specific energy £ min = — AE, returns. The time since disruption for this to happen is given 
by 



Ati = 2tiGM h (2AE)~ 3/2 m 0.068 yr f-^V" ( ^-Y (^) . (9) 



\WM Q J \M^J \R Q ; 

where we have used r P = r T and k = 3 2 . Material that is less bound takes a progressively 
longer time to return, and the overall light curve is determined by the distribution of debris 
mass as a function of binding energy. 

Rees (1988) assumed that the debris is uniformly distributed in mass between — AE 
and +AE. Numerical simulations (Evans & Kochanek 1989; Ayal, Livio, & Piran 2000) 
have shown that this is a reasonable approximation. The bound material then returns to 
pericenter at the rate (Phinney 1989; Evans & Kochanek 1989) 

2AM ft - t D \~ 5/3 A ft-t D ^ ~ 5/3 



M =3A^UfJ =A \^) AM = fM, , (10) 

where to is the time of the initial tidal disruption, AM is the actual mass that falls back to 
pericenter, which is a fraction / of the original mass of the star, and 

In Rees' original model, the whole star is disrupted and half the debris is on bound 
orbits, so that / = 0.5. However, a recent numerical simulation shows that not all the 
returning material is captured by the black hole: about 75% of the returning mass becomes 
unbound following the large compression it experiences on the way back (Ayal, Livio, & 
Piran 2000). This gives rise to a smaller / « 0.1. Another possibility giving rise to a small / 
is that the star is only partially disrupted: its envelope could be stripped by the black hole, 
leaving most of its core nearly intact (Renzini et al. 1995; see also §5 below). Here we treat 
/ as a free parameter. 



2 If the star is on a parabolic orbit with rp < tt, then the spin- up effect will be less important and we 
will have a smaller k (but k is still > 1), which will give rise to a At\ that is somewhat larger than that 
given by eq. (9). 
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The gravitational potential energy available from fallback is determined by the difference 
between the specific binding energy of the circularization orbit at r = 2rp = 2r? and the 
specific binding energy of the incoming material. Since Mh ^> M+, all the bound debris 
is on highly eccentric orbits with a specific binding energy much smaller than the binding 
energy of the final circular orbit. Thus, assuming that the fallback material radiates the 
energy release promptly, the radiation efficiency e is independent of time during fallback: 

GM+ fM H \ 2/3 ( M u \ 2/3 f MA 1/3 f , x 

The luminosity of the fallback process is then given by 

The luminosity peaks at t = t D + At 1; i.e. when the most bound debris falls back to the 
pericenter, so we have 

4. Comparison with the Flare in NGC 5905 

In §2 we have fitted the decline of the NGC 5905 light curve to a (t — t D )~ 5 / 3 law, and 
obtained the result given in equation (1) with a disruption time of t D = 1990.36 yr. The fit 
is acceptable (see Fig. 1), and the time at which the light curve reached its peak is fairly 
well constrained. Since the initial measurements at t ~ 1990.54 yr showed flux rising with 
time, the maximum in the light curve was clearly later than this time. At the same time, the 
maximum could not have been later than 1990.56 yr, since the fitted light curve in (1) falls 
below L 4.47 x 10 42 erg s _1 at that time. Thus, the light curve fit gives Ati within the 
range 0.18 to 0.20 yr, where At\ is the time lag between the peak of the light curve and the 
moment of tidal disruption. However, there is an error of ±0.03 yr on each limit. Therefore, 
we use conservative bounds on Ati 

0.15 yr < At x < 0.23 yr. (15) 

/,From the discussion in §2, the mass involved in the outburst of NGC 5905 is only 
~ 1O" 4 M . We therefore consider the following two possibilities for the candidate stellar 
object disrupted in NGC 5905: 
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(1) A low mass star (LMS), whose radius and mass are related by R± = R Q (M*/M Q ) 
(O.O8M < M* < 1M ). 

(2) A substellar objects (SSOs, including brown dwarfs and planets), whose radius and mass 
are related by R* = O.O6 J R (M^/M )- 1 / 8 (O.O1M < < O.O8M ) (Chabrier & Baraffe 
2000) 3 . 

Substituting these two radius-mass relations into equation (9) and noting the limits on 
At i in equation (15), we find that the black hole mass in NGC 5905 and the mass of the 
disrupted star are restricted to the two relatively narrow shaded bands shown in Figure 
2. The figure makes use of an additional constraint, namely a lower limit on the tidal 
radius tt- If the black hole is not rotating, then the smallest distance to which a star on a 
parabolic orbit can approach and yet not be captured by the black hole is r = 4r g , where 
r g = GMh/c 2 = th/2. If we have a maximally rotating black hole and if the angular 
momentum of the stellar orbit is parallel to the spin of the black hole, then the closest 
approach possible is r = r g . The real situation is probably in between these two limits, 
so we use a fiducial limiting radius of r = 2r g . Using equation (6), we have calculated 
relations connecting M H and M* for all three limiting radii and plotted them as dotted lines 
in Figure 2. 

The mass of the black hole in NGC 5905 is not very well constrained, but it is believed 
to be in the range M H ~ 10 7 — 1O 8 M (Komossa 2001), based on the correlation between 
bulge blue luminosity and black hole mass for spiral galaxies (Salucci et al. 2000) and on 
the correlation between bulge velocity dispersion and black hole mass for ellipticals and 
spirals (Gebhardt et al. 2000; Merritt & Ferrarese 2001). This information is sufficient to 
constrain the mass of the disrupted star fairly tightly: the mass must be either in the range 
0.01 - O.O2M in the case of a SSO or 0.6 - 1M in the case of a LMS. 

The result, however, depends sensitively on the choice we make for the tidal spin-up 
parameter k discussed in §2. Figure 2 has been calculated for k — 3, a reasonable and 
possibly likely value, but in principle k could be as small as 1. Figure 3 shows how the mass 
constraints change when we use the latter value. We see that the allowed range of models 
is limited to somewhat lower black hole masses, unless the black hole spins very rapidly and 
allows tidal disruptions down to r T ~ r g . We do not discuss the k — 1 case further. 



3 Chabrier & Baraffe (2000) state that their fit is good down to M± = 0.01M Q , but a comparison with their 
plot shows that the fit is acceptable even for lower masses. We apply their scaling down to M± ~ 0.006M Q . 
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^From equations (13) and (1) we infer that 

M H f MA 2/3 _ 0.0019 
1O'M {M e J ~ f ' Ubj 

where / is the fraction of the mass of the star that returns as fallback debris. The dashed 
lines in Figure 2 correspond to / = 0.0005, 0.012, and 0.12 respectively. We see that, for the 
case of a LMS, / must be < 0.0005, which means that an unusually small fraction of the 
star must have participated in the fallback. For the case of a SSO, somewhat larger values 
of / are obtained (but still rather small). 

In Figures 4 and 5 we plot the luminosities given by equation (13) for different values 
of M H) M*, and /, where the disrupted object is assumed to be a LMS (Fig. 4) and a SSO 
(Fig. 5), respectively, and the disruption time is taken to be t D = 1990.36 yr. We have 
chosen values of (M H , M+) from the shaded regions in Figure 2 so that equation (15) and the 
condition tt > 2r g are satisfied, and chosen / such that L pca k > 4.47 x 10 42 erg s _1 where 
-kpeak = L(t — to + Aii) is the peak luminosity. 

We confirm that, in order for the models to fit the observational data, a very small / 
(< 0.0003) is required for the case of a LMS (Fig. 4). For the case of a SSO (Fig. 5) a fairly 
small / (< 0.04) is again required if M H >2x 10 7 M Q . A largish / (> 0.1) can be obtained 
only if the black hole has a very small mass: M H < 4 x 1O 6 M . 

In Figures 4 and 5 we also indicate the ratio of the peak flux, F pea k, in the fitted light 
curve to the Eddington flux, i^dd, where -F pea k = L pca ^/ {-kR 2 x ) with Rx given by equation 
(4), and F Edd = L Edd / (4nr%,), with L Edd = 1.3 x 10 45 erg s' 1 (M H /10 7 M & ) and r T given by 
equation (6). We see that the ratio -F pea k/-^Edd is significantly larger than unity if the color 
temperature factor f c ~ 1, and it is not much less than unity even if f c ~ 3 (the likely value). 
We discuss the implications of this result in the next section. 



5. Summary and Discussion 

The X-ray light curve of the flare in NGC 5905 agrees with the (t — t_o)~ 5//3 dependence 
predicted for fallback of debris after tidal disruption of a star by a supermassive black hole 
(Fig. 1). If the black-body fit presented in Bade et al. (1996) is a good representation of the 
spectrum of the emission, then the bolometric correction to the X-ray flux is not very large 
and the observed light curve gives a meaningful measure of the energetics and luminosity of 
the flare. We note that, during and after the X-ray outburst, the optical brightness did not 
show any significant variation (Komossa & Bade 1999). This confirms that the bulk of the 
fallback emission probably did occur in X-rays. (If the flare had substantial hidden emission 
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in the far ultraviolet band, it is hard to see how there could have been no variations in the 
optical.) 

In addition to fallback, a second independent stage of evolution is expected when the 
circularized fallback material viscously accretes onto the black hole. This stage is expected 
to develop much more slowly and to be very much dimmer than the fallback stage (Appendix 
A). It is not relevant for the particular flare in NGC 5905 considered here. 

The time lag between the estimated disruption time t D = 1990.36 yr of the event and the 
time ti corresponding to the peak of the light curve gives the orbital period Ati of the most 
tightly bound debris. The data constrain Ati fairly tightly (eq. [15]). By the arguments 
described in §3, this translates into constraints on the mass Mh of the black hole and the 
mass M* of the disrupted star. The constraints are shown in Figure 2. 

In deriving equation (9), we have estimated the binding energy of the most tightly bound 
debris assuming that the star is fully spun up by tidal interaction before being disrupted. 
This corresponds to a spin-up parameter k = 3 (§3); this value is supported both by linear 
analysis (Press & Teukolsky 1977) and numerical simulations (Alexander & Kumar 2001). 
If on the other hand there is no tidal spin- up, then we expect k = 1, and the results change 
significantly (Fig. 3). We consider k = 1 to be somewhat unlikely, and have focused on the 
k = 3 results. 

There is also ambiguity in the value of the parameter in the definition of tt [eq. (5)]. 
Numerical simulations show that \i is usually of order unity, the value we use, but that it can 
be as large as 1.7 in specific cases (Frolov et al. 1994; Diener et al. 1997). The time interval 
Ati [eq. (9)] is very sensitive to the value of /i: Ati oc //. Because of the tight observational 
constraint on Ati (Figs- 2 and 3), fi cannot be much larger than 1 since otherwise it would 
imply too small a black hole mass. 

Another constraint on Mh and M* is obtained from the requirement that the radius tt 
at which the star is tidally disrupted must lie outside the radius of the marginally bound 
circular orbit. This requires r T > Ar g for a Schwarzschild black hole and r T > r g for the 
most favorably oriented orbit around a maximally spinning Kerr hole. The uncertainty in the 
limiting tt is indicated by the dotted lines in Figure 2. For concreteness, we have assumed 
that the limiting radius is intermediate between the two extreme cases mentioned above; we 
have used the limit Tt > 2r g . Yet another source of uncertainty derives from the fact that 
we have used Newtonian dynamics in our analysis, whereas it is clear that relativistic effects 
should be important (Luminet & Marck 1985; Laguna et al. 1993; Frolov et al. 1994; Diener 
et al. 1997). 

Ignoring these caveats, we obtain a tight correlation between the mass of the disrupted 
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star in NGC 5905 and the mass of the central black hole. As Figure 2 shows, the data are 
consistent with two distinct solution branches. 

One solution branch corresponds to a relatively massive black hole, M H close to 1O 8 M , 
tidally disrupting a main sequence star with M± ~ M . The other branch consists of 
a lower mass black hole, M H < 1O 7 - 5 M , disrupting a brown dwarf or a giant planet with 
M* < O.O2M . Assuming that NGC 5905 satisfies the well-known correlations between black 
hole mass and host galaxy properties (Salucci et al. 2000; Gebhardt et al. 2000; Merritt & 
Ferrarese 2001), the black hole mass is estimated to be M H ~ 10 7 - 1O 8 M (Komossa 2001). 
This mass range permits both solution branches. 

Figures 4 and 5 show sample theoretical light curves corresponding to the two solution 
branches. In each case, there is one free parameter, namely the fraction / of the mass of the 
disrupted star that returns as fallback debris to produce the observed radiation. We find that 
/ is unusually small if the disrupted star is a main sequence star; for the particular example 
shown in Figure 4, we require / < 0.0003, which is an extremely small value compared 
to the theoretically expected value of / ~ 0.1 — 0.5. Brown dwarf models give somewhat 
larger estimates for / (Fig. 5), though still uncomfortably small, unless the black hole has an 
unusually low mass (< 1O 6 - 5 M ). Though the estimated / can be increased by adopting a 
larger ambiguity factor /x in the definition of ty, we do not consider this likely since a larger 
/i will lead to an unacceptably small value of Mh- 

The above results on / are not surprising. As we showed in §2, the total radiative energy 
in the flare is only 4 x 10 49 erg, which for any reasonable radiative efficiency corresponds to 
an unexpectedly small mass ~ few x 1O~ 4 M . Figures 4 and 5 merely confirm this result 
using a detailed estimate of the radiative efficiency (eq. 12). 

Can / be so small for a tidally disrupted solar mass star? Rees (1988) and Renzini et 
al. (1995) have argued that a star could lose just its outer layers if the minimum distance 
from the star to the black hole is somewhat larger than tt- Here we point out a related 
possibility. The usual formula for the tidal disruption radius, tt = R*(Mh /M+) 1 / 3 , is correct 
only for a homogeneous star with a uniform mass density. For a real star, we could define a 
variable tidal disruption radius, 



r T (R) = R 



r m h - 


1/3 


" 3M H 


1/3 


. M (R)_ 




_47rp(i?)_ 





(17) 



where R < R* is a radius inside the star, M(R) is the mass contained within R, and p(R) is 
the corresponding mean density. The meaning of r T (R) is as follows: if the star is located 
at a distance r = r T {R) from the black hole, then the part of the star outside radius R 
(i.e., the spherical envelope between R and R+) is tidally stripped off, and the core inside 
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R remains intact. Obviously, r^max = rr(-R*) and rr, m in = (3iWjj /47rp c ) , where p c is 
the central density of the star. If the star is at a distance r from the black hole such that 
r T,min < r < 7Y,max, then only a part of the envelope of the star is peeled off. 

One important question is how the core of the star responds when it loses its envelope. 
What is relevant is the response on a dynamical time. Webbink (1985) has calculated the 
logarithmic derivative of the radius of a main sequence star with respect to mass, ( s = 
(din It*/ d\n M*) s , where the subscript s indicates that the entropy profile of the star is held 
fixed (as appropriate for a dynamical process). Consider a star that loses a fraction of its 
mass as a result of a tidal interaction and assume that it dynamically adjusts to a new 
equilibrium corresponding to its reduced mass. If ( s > 1/3, then the mean density of the 
new configuration will be higher than the original mean density of the star. The leftover 
mass is then stable to further stripping. However, if ( s < 1/3, the new configuration will be 
prone to further tidal mass loss, and one has a runaway situation in which the entire star is 
likely to be disrupted. 

Webbink (1985) shows that ( s > 1/3 for \og(MjM Q ) > -0.2. Since this is precisely the 
mass range of interest if the disrupted star in NGC 5905 was a main sequence star (Fig. 2), 
we conclude that partial stripping is a viable possibility. To estimate what fraction of its 
mass the star might lose, consider a star like the Sun as an example. For the Sun, the mass 
inside a radius R± = 0.7R & is 0.975M Q , and correspondingly, tt{Ri) = 0.706 rY jmax . Thus, 
for a range of distances between rj> imax and 0.77Y max (a non-negligible range), such a star 
will lose only a couple of percent of its mass. If we assume that 10% of this mass returns via 
fallback, then we can have / as small as ~ 0.002. This is still much larger than the value 
needed to explain the flare in NGC 5905, but it is at least within striking distance. 

Fully convective lower main sequence stars and degenerate brown dwarfs have ( s ~ —1/3. 
These stars are likely to be fully disrupted in tidal encounters since they expand to a larger 
radius as they lose mass. 

Apart from the small value of /, another unexpected result of our analysis is that the 
peak radiative flux emitted during the flare, -F pca k = o~T^ s = a(T^/ f c ) 4 , is close to the 
Eddington flux. (Here T bb is the black-body temperature of the radiation and f c is a color 
correction factor.) The results are shown for our models in Figures 4 and 5. Even with 
fc ~ 3 (§2), we find -F pea k/-pEdd ~ 0.1 —0.3; the value increases rapidly with decreasing f c (as 
fc 4 )- What makes this surprising is that the luminosity of the source, even at the peak of 
the flare, is much below the Eddington luminosity: compare L peak ~ 4.5 x 10 42 ergs -1 with 
-^Edd ~ 10 45 (M#/10 7 M Q ) ergs -1 . Clearly, the radiation in the flare must have been emitted 
from a small patch on the source, and the patch must have projected a small solid angle at 
the central black hole. This result was anticipated in the model-independent analysis of §2, 
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where we showed that the effective size of the emitting patch is small compared to the area 
we might expect if the entire accretion disk out to a few r g is involved in the emission. 

The small size of the emitting region can be understood in analogy with mass transfer 
binary stars. In the latter systems, it is known that the cold gas stream from the secondary 
star impacts the accretion disk at a compact "bright spot" on the outer edge of the disk. 
Most of the radial kinetic energy of the stream is converted to thermal energy near this 
spot and radiated locally. The analogy with the tidal disruption problem is close because 
the returning material after a disruption event comes in similarly as a narrow cold stream 
(Kochanek 1994; Kim, Park, & Lee 1999). It is therefore natural that the radiation during 
the fallback stage should be emitted from a relatively small patch on the source. The effect 
is probably even more severe in the tidal disruption problem. Lense-Thirring precession may 
well cause the incoming stream to be misaligned with respect to the orbital plane of the 
accretion disk (Kochanek 1994), causing the geometric area of the impact spot to become 
smaller. (We note that in binaries the stream appears to skim over the surface of the disk 
after impacting at the bright spot, thereby increasing the effective area over which the kinetic 
energy of the stream is dissipated. This effect would be reduced if the orbit of the stream is 
not in the plane of the disk.) 

Having decided that emission from a small patch is natural for this problem, we should 
next ask whether it is a pure accident that the flux from the patch in the NGC 5905 flare 
happens to be nearly equal to the Eddington flux. We would like to suggest that it is perhaps 
not an accident. 

Imagine a narrow stream with a mass accretion rate M impacting on the accretion torus 
at a radius tt and over a small patch that subtends a solid angle Af2 strC am at the black hole. 
Let us suppose that the kinetic energy in the stream is thermalized and radiated immediately 
from the patch. The radiated flux is locally given by eMc 2 /Af2 strcam r|,. What would happen 
if this flux is greater than the local -FEdd? Jeremy Goodman (private communication) suggests 
that the super-Eddington flux would cause the local radiating atmosphere to puff up and 
that the radiation would be emitted from a larger patch than the original impact region. 
Moreover, it would be natural for the radiating patch to self-adjust so that the flux is close to 
the Eddington limit. According to this picture then, the radius Rx estimated in equation (4) 
refers to the size of the puffed up atmosphere, while the actual transverse size of the stream 
is smaller. If the stream is extremely narrow, then the flux would be close to Eddington 
during a significant fraction of the light curve. Equivalently, the black-body temperature of 
the radiation would remain constant, even while the luminosity decreases. 

This picture (which is predicated on the assumption of a narrow fallback stream) 
suggests that most tidal disruption events should be visible in soft X-rays; the condition 
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F = F E dd, coupled with f c ~ 3 (see §2), ensures that the blackbody temperature of the radi- 
ation will be > 50 eV for a 1O 7 M black hole, which is well-matched to X-ray observations. 
Soft X-ray surveys should thus be well-suited for an investigation of both the frequency of 
tidal disruption events and the amount of mass accreted in each event. In this connection, it 
is interesting that, with the exception of the UV flare of Renzini et al. (1995), all known can- 
didate flares for tidal disruption events have been identified in the X-ray band (see Komossa 
2001 for five candidates and Choi et al. 2002 for a possible sixth candidate). 

Menou & Quataert (2001) have emphasized that, given current estimates of the rate 
of tidal disruptions, the nuclei of many nearby normal galaxies should be much brighter 
than observed if a fair fraction of a stellar mass is accreted by the black hole after each 
tidal disruption event. This is not relevant for the late accretion phase of the particular 
NGC 5905 flare discussed here because, as we have discussed (see also the Appendix below), 
very little mass (<C 1M Q ) returned to the BH in this case. However, whether we attribute 
the small amount of mass in the NGC 5905 flare to partial stripping of a star or to the 
disruption of a substellar object, it does not help to explain the general problem of why 
nearby galactic nuclei are so faint. In addition to partial disruptions of the sort that might 
have occurred in NGC 5905, we expect total stellar disruptions in significant numbers and 
these should provide a significant mean mass accretion rate in galactic nuclei. Disruptions 
of substellar objects only add to the purely stellar disruptions considered to date. Thus, 
alternative ways of explaining the dimness of nearby galactic nuclei are required (Menou & 
Quataert 2001; Narayan 2002), beyond the specific scenario outlined here for the giant flare 
of NGC 5905. Either the mass must accrete in a radiatively inefficient mode or the accretion 
must be suppressed or slowed down significantly. 
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A. Viscous Evolution of the Thin Disk Formed from Tidal Disruption 

After circularization, the fallback material from a tidal disruption event will form a 
torus at a radius r ~ 2rr- The torus will then spread out viscously and accrete onto the 
black hole. 

The time-scale for this stage of the evolution is determined by the viscous time-scale 
(Cannizzo, Lee, & Goodman 1990; Ulmer 1999) 

t ™ = J ^ K2 - 9xl0yr (al) \W*) U) U) • (M) 

where txe P is the Keplerian orbital period, a is the standard viscosity parameter (Shakura & 
Sunyaev 1973), and h is the ratio of the disk height to radius. It is easy to check that for the 
case of the flare in NGC 5905, where L pcak ~ lO~ 3 LEdd(^fe/lO 7 M ) _1 , the accretion disk 
would correspond to the regime of the "middle disk" in which electron scattering opacity 
dominates over the Kramers opacity and gas pressure dominates over radiation pressure 
(Shakura & Sunyaev 1973; Novikov & Thorne 1973; Frank, King & Raine 1992). Then, h is 
estimated to be 

h „ lx 10 -3 (JL)-*" pM 1 " (4M" 1/10 f-f ° • (^) 

\0.lJ \lO-*M md ) \10 7 M & ) \r H J 1 J 

where M 2 is the mass accretion rate, Msdd = -^Edd/(0.1c 2 ), and r is the radius in the disk. 
In equation (A2) we have neglected a dimensionless factor ~ I. 

Inserting equation (A2) into equation (Al) and setting r = 2r T , we obtain 

-2/5 

t vis w 2.3xl0 3 yr 

Ki^) \M) ■ (A3) 

The large value of t vis confirms that the flare in NGC 5905 corresponded to the fallback 
stage, not the viscous accretion stage. 

From the discussion in §3 on the fallback and viscous accretion stages associated with a 
tidal disruption event, the ratio of the total energy radiated during the accretion stage to the 
total energy radiated during the fallback stage is 2r P /r ms — 1. Therefore, we can estimate 
the ratio of the peak of the viscous luminosity to the peak of the fallback luminosity by 

£yis _ ( 2r P ^ At! 
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V0.17 VlO- 3 ; \WM G ) \M Q J \r ms J y ' 
where equations (9) and (Al) have been used. 

The extremely large value of t v ; s and the extremely small value of L vis /Lft, imply that 
the accretion stage will never be observed in practice for the particular case of NGC 5905. 
[In fact, according the above estimates, one may expect a new tidal disruption before the 
disk evolves significantly (Syer & Ulmer 1999; Magorrian & Tremaine 1999).] 
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Table 1. ROSAT Data on the NGC 5905 Flare 



TIME a CR b ERR_CR C LUMX d ERR_LUMX' 



1990. 


53 


2.2981E- 


1 


8.1771E- 


-2 


1.72 




0.61 


1990. 


.54 


2.8713E- 


1 


7.3284E- 


-2 


2.15 




0.55 


1990. 


.54 


4.1855E- 


1 


6.2263E- 


-2 


3.14 




0.47 


1990. 


.54 


4.5519E- 


1 


7.5933E- 


-2 


3.41 




0.57 


1990. 


.54 


4.1786E- 


1 


6.2403E- 


-2 


3.13 




0.47 


1990. 


.54 


5.9612E- 


1 


8.7399E- 


-2 


4.47 




0.66 


1990. 


.54 


5.4012E- 


1 


8.2004E- 


-2 


4.05 




0.62 


1990. 


.97 


<0.09 








<0.68 






1991. 


.04 


<0.10 








<0.75 






1993. 


.55 


7.0000E- 


3 


1.0000E- 


-3 


4.00E- 


-2 


5.71E 


1996. 


.89 


2.9000E- 


3 


4.0000E- 


-4 


1.66E- 


-2 


2.29E 



a Time of measurement in years. 

b Count rate in cts s _1 . 

c Error in the count rate. 

d X-ray luminosity in units of 10 42 erg s _1 . 

e Error in the X-ray luminosity. 

Note. - - The data were provided by S. Komossa. The con- 
version from count rate to luminosity was done assuming a 
black-body spectrum and a Hubble constant of H = 50 km s _1 
Mpc -1 . For this model, 1 cts s _1 corresponds to ~ 7.5 x 10 42 
erg s _1 during 1990 when the source was in the high-state, and 
~ 5.7 x 10 42 erg s _1 during 1993 and 1996 when the source had 
a flatter spectrum (S. Komossa 2001, private communication). 



- 21 - 




Fig. 1. — The X-ray light curve of the flare in NGC 5905 (see Table 1 for the data), fitted to 
a model of the form (t — t D )~ 5 / 3 , as expected for fallback of material after a tidal disruption. 
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Fig. 2. — Restrictions on the mass of the black hole, M H , and the mass of the disrupted 
stellar object, M*, derived from the X-ray flare of NGC 5905. The shaded areas show the 
allowed regions. It has been assumed that the star is spun up to breakup at pericenter 
(k — 3). The region to the left of the thin vertical line corresponds to brown dwarfs/planets, 
and the region to the right to low mass stars. The limits used are (1) 0.15yr < At x < 0.23 yr, 
where Ati is the time for the most bound disrupted material to fall back to pericenter; and by 
(2) rx > 2r g , where tt is the tidal disruption radius (assumed to be equal to the pericentric 
radius) and r g = GMh/c 2 . For comparison, lines corresponding to tt = r g and tt = 4r s are 
also shown. The dashed lines correspond to equation (16) for various values of the fraction 
/ of the mass of the disrupted star that falls back. 
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Fig. 3. — Same as Fig. 2, but assuming that the star experiences negligible spin-up (k 
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Fallback model: t D = 1990.36, 
M H = 8 xlO 7 M Q , M + = 1 M Q , 
(f. F peak /F Edd ) = (0.0003, 0.34 f-J), 
(0.00026, 0.30 f^ 4 ), 
(0.00023, 0.26 f" 4 ), 
(0.0002, 0.23 f" 4 ) 
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Fig. 4. — Comparison of light curves predicted by the fallback model with data: the disrupted 
object is assumed to be a low mass star. The parameter / is the ratio of the accreted mass 
to the total mass of the star, F pe3 ^/ F^d is the ratio of the peak flux = L pca k/ \-kB? x ) to the 
Eddington flux = LEdd/(47rrf,), where / C 3 = / c /3. Four models are shown, each with the 
values of / and Fpeak/^Edd indicated (from top to bottom). 
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Fallback model: t D = 1990.36, 



M H = 2 xlO 7 M , M„= 0.016 M , 

C Fpea k / F E d d) = (0-04, 0.19 f-J), 

(0.032, 0.15 f" 3 4 ), 
(0.025, 0.12 f" 3 4 ), 
(0.019, 0.09 f- 3 4 ) 




Fall-back model: t D = 1990.36, 
M H = 3.2 xlO 6 M , M + = 0.0075 M G , 

(f. Fpeax/FEdd) = (0-2, 0.34 f c - 4 ), 
(0.18, 0.31 f- 4 ), 
(0.16, 0.27 f- 4 ), 
(0.14, 0.24 f- 3 4 ) 
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Fig. 5. — Similar to Fig. 4, but assuming that the disrupted object is a brown dwarf or a 
planet. Two cases are shown: Mh = 2 x 10 7 M q , brown dwarf of mass M* = 0.016M Q (upper 
panel); and M H = 3.2 x 1O 6 M , planet of mass M* = O.OO75M (lower panel). 



